Second order gravitational effects on CMB temperature anisotropy 

in A dominated flat universes 



Kenji Tomita 

Yukawa Institute for Theoretical Physics, Kyoto University, Kyoto 606-8502, Japan 

Kaiki Taro Inoue 

Department of Science and Engineering, Kmki Universtty, Higashi-Osaka, 577-8502, Japan 

(Dated: May 15, 2008) 

We study second order gravitational effects of local inhomogeneities on the cosmic microwave 
background radiation in flat universes with matter and a cosmological constant A. We find that 
QQ ' the general relativistic correction to the Newtonian approximation is negligible at second order 

' provided that the size of the inhomogeneous region is sufficiently smaller than the horizon scale. 

, For a spherically symmetric top-hat type quasi-linear perturbation, the first order temperature 

^Nj ■ fluctuation corresponding to the linear integrated Sachs- Wolfe (ISW) effect is enhanced(suppressed) 

by the second order one for a compensated void(lump). As a function of redshift of the local 
inhomogeneity, the second order temperature fluctuations due to evolution of the gravitational 
potential have a peak before the matter-A equality epoch for a fixed comoving size and a density 
contrast. The second order gravitational effects from local quasi-linear inhomogeneities at a redshift 
2; ~ 1 may significantly affect the cosmic microwave background. 
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PACS numbers: 98.80.-k, 98.70.Vc, 04.25. Nx 



dl' I. INTRODUCTION 

O 

' Generation of temperature anisotropy in the Cosmic Microwave Background (CMB) due to linear pertubatons 
of gravitational potentials is called the Sachs- Wolfe (SW) effect^. There are two contribution to the SW effect. 
The first one called the non-integrated Sachs- Wolfe effect is produced by fluctuations of gravitational potentials 
at the surface of last scattering. The second one called the integrated Sachs- Wolfe (ISW) effect is generated by 
04 ■ time-varying gravitational potentials as the CMB photons pass through them from the recombination epoch to 
^ [ the present epoch. 

• Recently, much attention has been paid to the non-linear version of the ISW effect called the Rees-Sciama 
■ (RS) effect [3] on cosmological scales, since the RS effect of local inhomogeneities at quasi- linear regime may 
. explain the origin of observed anomalies d, 0, [E B 0, Si in the large-angle CMB anisotriopy [§, [l3| . It is argued 
• ■ that the anisotropy for compensated asymptotically expanding local voids can be larger because the second 
, order effect enhances the ISW effect [11[. For compensated local lumps, the nature of the second order effect 

^ ■ at the A dominated epoch remains unknown. 

[ So far, two types of treatment have been used for deriving the RS effect of local inhomogeneities. One is the 
• • ■ general relativistic treatment [l^, [3, [3 , in which a non- linear version of the SW effect at second order has 
been studied. Another one is the simplified Newtonian treatment [l^, [l^, [13, EM [3 HO) [Ml in which influences 
of nonlinear halo clustering upon the CMB temperature fluctuations and gravitational lensing phenomena have 
been explored. 

In this paper, we first study the RS effect comparatively using the general relativistic second order pertur- 
bation theory and the Newtonian approximation and we show the consistency between the two approaches in 
the perturbative regime. Next, we consider the RS effect of a compensating quasi-linear void/lump modelled 
by a spherically symmetric top-hat type density perturbation at a redshift z and we investigate the temporal 
change in the linear and second order temperature fluctuations as functions of z. In §2, we derive the solutions 
of Einstein equations for general relativistic perturbations at second order in the Poisson gauge (generalized 
longitudinal gauge) and consider the limit of k (~ the scale of inhomogeneities / Hubble radius ) ^ 1. In §3, we 
study the perturbations in cosmological Newtonian approximations and their relation to the general relativistic 
perturbations up to second order. In §4, we study the temperature anisotropy owing to first and second order 
perturbations based on relativistic perturbation theory 22]. In §5, we investigate the correlations between 
first order and second order temperature fluctuations and their temporal behavior for a spherical top-hat type 
density perturbation. §6 is dedicated to concluding remarks. In Appendices A and B, the main components of 
Einstein equations and the derivation of their solutions are shown. In Appendix C, the integrations of metric 
perturbations for spherically symmetric top-hat type perturbations along the light paths are shown. 
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II. GENERAL RELATIVISTIC SECOND ORDER PERTURBATIONS 



In what follows, we use the units of SttG = c = 1, the Greek and Latin letters denote 0, 1,2,3 and 1,2,3, 
respectively. Index "B" represents the value for the unperturbed background spacetime. Sij{= = S^^) are the 
Kronecker delta, and subscripts (n) correspond to n-th order quantities. 

As a function of conformal time fj and Cartesian coordinates £*, the metric of spatially flat Friedmann- 
Robertson- Walker (FRW) universes with first order and second order perturbations ip^^''\ 4'^"\ Xi]^\ n — 1,2 
is described by 



- (1 + 2?/^^^^ + ilP'^)d7f + 
(l-20W-0(^))5,+x!j^ + ix!f 



.(1) 



dx'dx^ 



(2) 



}■ 



dfjdx^ 



(2.1) 



where a{fi) is the scale factor, 77 (= x^) is related to the cosmic time t hy dt ~ a{fi)dfi, and Xif' satisfy 

= ^^"^Xhl ~ 0- In what follows, we use the normalized scale factor a — ciHq^, the comoving coordinate 
X* = x^/H^^, and the conformal time rj = fi/H^^ where H^^ is the Hubble radius at present 770. The scale 
factors at present are defined as ao = 1 and ao = . Note that the scale fatcor a has a dimension of length 
while 77 and r do not have dimensions in our notation. 

Energy density and 4-velocity of dust matter are written in terms of the background quantities and pertur- 
bations as 



(2.2) 



1 



1 



,(2)m 



where denotes the 3- velocity of the dust matter. From the condition g^iyM^'u'^ 



(2.3) 

-1, we obtain the relations 
(2.4) 



«(2)o = _^(2) + 3(^(i))2 ^ [2zf ) + z;W]7;(i)\ (2.5) 
Einstein equations for the unperturbed background with dust matter and a cosmolological constant A are 



iia' /of = {pB + PA)a^ 



(2.6) 



6(a7a)' = -(pB-2pA)a2 



(2.7) 



where pA is the energy density of the cosmological constant A and a prime denotes derivative with respect to 77. 
To fix the gauge freedom of perturbations, we adopt the Poisson gauge defined by 



5'™4"i = and <5'"xi:i = 



for n — 1 and 2. This gauge is a generalized version of the longitudinal gauge which is defined by zf^^ = 
and Xij^ = [23| . and gives us a metric expression convenient for a Newtonian interpretation, as well as the 
longitudinal gauge [2^ . 

In what follows, we consider only scalar-type perturabtions at linear order. Then, we have 



zW=0 and xil^=0 



(2.9) 



The Ricci tensor, the Einstein tensor and the energy-momentum tensor for dust matter are shown in Appendix 
A. Solving the Einstein equation we obtain the expression of first order scalar-type perturbations in the growing 
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mode in terms of functions P{ri) and f (x) as"'^ 



(2.10) 



(2.11) 



Sp^'^/PB = -^[ia'/a)P' - l]AF+^{a'/a)P'F, 



(2.12) 



-^[{a'/a)P'-l]F, ^(D^^ip'F. 



(2.13) 



where is dF/dx\ A is the Laplacian d / dx^ dx\ and 



V = I daa ^^^[n„M + nAQa^] 
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(2.14) 



where rimo and f^Ao are the density parameters for matter and the cosmological constant at present. P{'q) is 
the solution of the growing mode in equation 



P" + ^p' -1 = 0. 



(2.15) 



Note that the potential function P(x) is related to the first order matter density contrast £,„ in comoving slices 
(defined in the comoving synchronous gauge) as 



Psa" 



:[{a'/a)P' ~l]AF 



(2.16) 



Next we consider relativistic second order perturbations corresponding to the first order perturbations in the 
growing mode. From Eqs. (jB9p - (jB12() in Appendix B, which are derived by solving the Einstein equations 
52 =62^, we have 



b^^^ = CiFiFj + (2 ■ 100*0 + CsF^ + Ci ■ lOOOo, 



(2.17) 



V'^'^ = ^iFaF^i + 6 ■ 100*0 + CsF^ + U ■ lOOOo, 



(2.18) 



where 



Ci 

C3 



ip(l-^p' 

4 V a 
4 I a 



a 
a 



(2.19) 



— Ci, C2 — C2, 

6 = J{4-7^P' + 
4 L a 



,11 



a /a 
- + 5(- 

a \ a 



'\ 2 



u^'-h-'-^p' 

^ 6 1 a 



2a" /"^2 

+8 

a \ a 



(^)'](P')^},(2.20) 



^ In paper [12| | (referred as Paper [T3|). the first order perturbations in the Poisson gauge were derived by transforming the solution 
in the comoving synchronous gauge to that in the Poisson gauge [23I. |25|| . After the publication several misprints were found in 
Eqs. (4.6) - (4.8) of Paper [T3|. which should be taken into account for deriving expressions l|2.10ll - l|2.13| l. Perturbations in the 
decaying mode were derived in Paper [T3l . but they are omitted here. 
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.(2) 



P'(1 + P")C,, 



(2.21) 



and 



where C,- satisfies 



(2) 
V = 



200 



-iFiFi),,-F,AF 



and (5xy in Eq. p.22p satisfies 
^9772 

^1^0 and 8o are defined as 



2a' d 
a drj 
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" 200 



-A)5x^, = ^P'A'D,, 



1 r 
7 



p_^(p')2 



[AFf 



and (3(77) satisfies 



2a' 



Q' 



AOo = ^-0 



2^ ' . 



100 



(2.22) 



(2.23) 



(2.24) 



(2.25) 



(2.26) 



The above second order solutions are consistent with those shown in Eqs.(4.12) - (4.15) in Paper[T3|, which was 
derived using a transformation from the comoving synchronous gauge to the Poisson gauge. Note that ADij 
and A^D,, correspond to and in Eq.(2.25) of Paperfl^, as AAj = -G{ and A2A,, = -G{. Eq.(2.17) 
of Paper JJj with misprints must be replaced by the above correct equation (|2.26p . 

In the above solutions, the ratios of terms including 80 to terms including PjP; and ^'o are of the order of 
At^, where k = |x|/ry, and |x| and rj are the characteristic spatial scale of local perturbation and the horizon 
size. Therefore, the terms including are negligible for local perturbations that are sufficiently smaller than 
the horizon size. Thus for k ^ 1, we have 



0(2) = ijj 



(2) 



Ci PiPi + C2 100*0, 



(2.27) 



zf^ -P'(1 + P")C„ 



and 



(2.28) 
(2.29) 



III. COSMOLOGICAL NEWTONIAN APPROXIMATION 



In the cosmological Newtonian approximation, we assume that e = a|v|/c ^ 1 and k = |x|/7? ^ 1 but 
X = {p/pB ~ 1)^/2 is arbitrary, and consider only ip and 4>(= ^ the metric perturbations [2q . Then, from 
the difference between the perturbed equation Pg = Pq — and the background counterpart, we obtain 



AV' = ]^a^{p-pB)- 



(3.1) 



From the conservation equation T.^" = and the energy-momentum tensor (for the perfect fluid with pressure 
p ) Tf^" — {p -\-p)u^u'^ +pg^^'^, we obtain the equation of continuity 



and the equation of motion 



; 3a' , 
P + — p+ pw* = 
a 



(3.2) 



(3.3) 
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where = au^. By solving these equations (|3.ip - (|3.3p . we have p,v'^ and -0, which can determine the lowest- 
order RS effect in the form of spatial integration of V'.i along a light path. 

In what follows, we only consider the perturbative case with e ^ 1, k ^ 1, and x ^ 1- First, expressing the 
perturbations as 

i, = ^W + p = PB+5p^'^ + ^Sp<'^\ =v(^r + lv^^)\ (3.4) 

the first order equations for these perturbations in the pressureless case are given by 

Spi^y + ^5p('^+pBv'^'^\ = 0, 
a 

v'^^y' + ^^WVv'/'=o, (3.5) 
a 



and the first order solutions are 



Sp^^^/PB = ^[{a'la)P' ~ 1]AF, (3.6) 
= ip'F,, (3.7) 



^^'^=-\{\--P'\f^ (3.8) 



and 

2 V a 

where the above Sp^^^ , w*^^^' and '^'^^^ are equal to Eqs. (|2.12p . (j2.13|l . and (|2.10|1 in the limit of k ^ 1, respectively. 
Next the corresponding Newtonian second order equations are 

a 



+ 2 + %(2)V^(f)=0. (3.9) 



By substituting the above first order solutions Eqs. (|3.6|) . (|3.7p . and ()3.8|) to Eq. (|3.9p and eHminating ip''^^ and 
v^'^y , we obtain 

(Sp^'^/pb)" + -{Sp^'^/pb)' Ja^ V = \{P'f^{FiFi) - ^^[a'p'(^P' l)]\/\FF,),. (3.10) 
a 2 4 psa a 

Assuming that V'*-^' is given by Eq. (|2.27p and using the first line of Eg. p. 91) . we can express Sp'^'^^ as 

V/PB = ^{CiA(F,i^,) + |C2[F.,F., - (AF)2]}. (3.11) 
Pb^ z 

By substituting this Sp'^'^^ /pB to Eq. p.lOp . we find that it is a solution of Eq. (|3.10p . so that i/;'^' given by 
Eq. (|2.27p is also a solution of Eq. (|3.9p . Thus we proved that in the cosmological Newtonian limit with e <C 1, 
K <C 1, and X <C 1, the general relativistic second order solution is consistent with the Newtonian one. 

IV. TEMPERATURE ANISOTROPY 

In this section we consider the observed temperature of the CMB radiation which was emitted at the recom- 
bination epoch and received at the present epoch. The relation between the emitted and received temperatures 
Te To is 



To = (Wo/We)Te, 



(4.1) 
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FIG. 1: Time dependent terms and —("2 in second order temperature fluctuations as a function of redsliift z. f2mo 
denotes tlie present value for the matter density parameter. 
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where uj = —g^i,u'^k'^, m'' is the observer's and emitter's velocities, and fc^ (= dx^ /dX) is the wave vector of 
photons with affine parameter A, which satisfies the null geodesic equation in the perturbed FRW universe. 
Solving this equation, the first order and second order perturbations of observed temperature AT^^^ and AT^^^ 
were derived in the gauge-invariant manner by MoUerach and Materresefl^. When the background null geodesic 
rays are expressed as x^'^^'^ = (A, (Aq — A)e*) and k^^^^^ — (1, — e'), the first order temperature fluctuation is 



he 



(4.2) 



where 



he 



dXA^^y, 



(4.3) 



is the unit (three-dimensional) directional vector, the subscripts e and o denote the epochs of emission and 
observation, r is the temperature fluctuation at the emission epoch, and it is assumed that x^^^'^^Xo) = and 
k'^^^''{Xo) = 0. The integral term lie represents the contribution due to the ISW effect. 

Similarly, second order temperature fluctuation is given by 



AT^^^ /T 

he 



he + [he? + (AT(2)/T), 



(4.4) 



where (AT'-^)T/T)oe is the sum of terms consisting of second order and products of first order quantities at 
observer's and emitter's positions, and 



^(2) = ^,(2) + 0(2) + 



(2) 
z, e 



(4.5) 



For quasi-linear perturbations with £,„ = 0(0. 1), the contribution due to the RS effect can be written as 

f2) ('2') 

(AT/T)]^g — lie + he- For scalar-type hnear or quasi-linear perturbations, the terms ' and xlj in Ea. (|4.5p 
can be neglected. Therefore, the second order contribution l2e to the RS effect can be evaluated using only 
perturbations (j)^'^^ and ^'■-^"^ which are obtained from Eq. (j2.27p with the following conformal time derivatives 



4 L a 



(4.6) 



C2' 



187 



-P 



-P' 



1 

21 



pp' 



5 a' 
36 a 



{P- 



J\2 



(4.7) 
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FIG. 2: A light path in the spherical model. 
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FIG. 3: The matter density contrast for a top-hat type spherical void. 
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(2) (2) 

where and xlj are neglected. For perturbations sufficiently smaller than the Hubble scale, i.e., k = |x|/?7 ^ 

1, the dominant terms in and t/i^^^' are those multiplied by FjF^i and ^Pq. 

For the Einstein-de Sitter (EdS) model (A = 0), we have a ocif ,P = if /IQ and Q = 0, so that 

Ci = — C2 = ^r;2, and Ci'/Cz = - — ■ (4-8) 

200 ' ' 210 ' ' ^^'^^ 20 ^ ' 

For a A-dominated model with (rimOj ^^Ao = (0.27, 0.73), numerical calculations of Ci and C2 show that C1/C2 have 
the values -2.68, -2.99, -3.08, -3.12, -3.14, -3.15 for the redshifts z = 0.05, 0.1, 0.2, 0.5, 1.0, 2.0, respectively. 
For z ^ 1, C1/C2 is nearly equal to —63/20 in the EdS model. As shown in Fig. [U C( and C2 have a peak well 
before the matter-A equality epoch z„ja = ((1 — ri,„o)/^^mo)^^'^ — 1, implying that the second order contribution 
to the RS effect in the quasi-linear regime is not so important at an accelerating epoch. 



V. BEHAVIORS OF FIRST ORDER AND SECOND ORDER TEMPERATURE FLUCTUATIONS 

IN A SPHERICAL TOP-HAT MODEL 

In order to investigate the nature of the temperature fluctuations exerted by local inhomogeneities, we consider 
a simple toy model with a spherically symmetric density perturbation. In what follows, we assume that we are 
outside the local density perturbation (Fig[2]). In this model, the potential function F(x) can be written as a 
function of the comoving distance r = [x^]^/^ from the center of the perturbation as F = F(r). Then, we have 
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(2) (2) 

Ci = 0, Dij = and = xlj = 0- For the top-hat type matter density perturbations, the functional form for 
F{r) is given in terms of constant parameters b and c as 

^P-^i(^'Fr)-c, -h, (5.1) 

for < r < 7'Oj To < f 1^ fi^ fi < r, respectively(Fig. [3]). This model represents a void if c > and 6 > 0, or a 

lump if c < and 6 < 0. Moreover, if cro^ = b{ri^ — ro^), the mass is totally compensating. —AF represents 

a value that is proportional to the matter density contrast in a comoving gauge. Note that —F{x) describes 

the X dependence of the gravitational potential ip. Here, tq and ri are the inner radius and the outer radius, 
respectively, ri — tq corresponds to the width of the wall. 
Integrating Eq. (|5.ip . we obtain 

1 r 

Fr^ ^ drr^AF. (5.2) 
Jo 

Then, under the condition that F,r is regular at the center r = 0, we have 

P,r = lcr, l{c + b)riy-^br, ^[c^ ^ birf - rl)] (5.3) 
for < r < ro, ro < r < ri,ri < r, respectively. From Eq. (j2.25|) . on the other hand, we have 

-i^(r2vE'o,.),. = [r(F,)2],„ (5.4) 

1 DO f^(F)^ 
- —IMr)) - MO)] = ^-^dr ^ L (5.5) 

Under the condition that ^'o is regular at the center, ^'o(O) is determined from the boundary condition that the 
perturbation is local, i.e., ^'o — ^ as r — > oo. The explicit forms of / and ^'o(O) are shown in Appendix C. 
Integrating Ea. (|5.3p . we obtain 

F = --C (r^ - rh +Fo for < r < tq 
6 

= --(c + 6) r3(- - — ) - -6 (r^ -r2) + Fo for tq < r < n 

^ [crl - b{rl - rl)] for ri < r, (5.6) 



so that 



3r 



where 



Fo 



\{c + b) rlQ- -]r)+\b (rf -rl)- " Hrl - rl)]. (5.7) 



As shown in Fig. |31 the thinner wall (i.e.,(ri — ro)/ri ^ 1 ) gives deeper potential for a lump provided that the 
matter density at the center is fixed (i.e., c is constant). 

Using Eqs. (|2.10p . (|4.2p . and (|4.3p . the first order temperature fluctuation due to the local perturbation is 
given by the integral term in Eq. (|4.2p . which represents the ISW effect, 



(AT(i)/ 

where 









\e — 


-2 










-a" 


-K 


a'\2- 


- a 




a / . 



(5.8) 



Now we consider a light path passing through the center of a spherical compensating void/lump and assume 
that its outer radius ri is sufficiently small compared with the Hubble radius at the time the photons pass 
through the center. Then the integration of F derived in Appendix C reduces to 



Fdr = --c(ro)3 (1 + y) ln(l + l/y) 



(5.10) 
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where y = b/c and ri/rg = (1 + \/y)^^^. From the above equations, we can express the first order temperature 
fluctuation as 



-c{T^fw^{y) - + _-3 - P' 
9 La \ a \ a / / 



n=Va 



(5.11) 



where wi {y) — —y ln(l + 1/y) wliich is a negative definite function for y > and rjQ — rjc represents the comoving 
distance to the void/lump. Using Eqs. (|2.16p . (|2.6p . and (|2.7[) . c can be written in terms of the matter density 
contrast {em)c in the comoving gauge at the center of the inhoniogeneity as 



(e„), 2[2(^) 



a 

a 



-P' -I 



(5.12) 



From Eq. (|5.1ip and Ea. (|5.12p . we find that the first order temperature fluctuation due to a compensated void 



or lump is approximately written as 0[{/S.T^^^ /T) 



loc\ 



(^m)c 



As shown in Figs. [5] and [H] a compensated 



spherical void redshifts the photons whereas a compensated spherical lump blueshifts the photons at first order 
irrespective of the width of the wall. 

In a similar manner, we can express the second order temperature fluctuation when the photons pass through 
the center of a compensating void/lump {— ~\ fx" dXA^'^^ ) as 



(5.13) 



'W2{y) = y[l - 2/ln(l + l/y)], 



(5.14) 



where c is given by Eq. (|5.12p . Note that W2{y) is a positive definite function for y > 0. The order of the second 
order temperature fluctuation is OKAT^^^^ /T)ioc] ~~(em)c^i- As shown in Figs. [5]and[6l either type (void or 
lump) of density perturbation redshifts the photons irrespective of the width of the wall. This behavior suggests 
that the second order gravitational effect leads to a flow of matter from the wall to inside the wall in either case. 
Then the gravitational potential becomes smaller and photons passing through the center of the void(lump) get 
further redshifts. 

To see this quantitatively, we calculate the ratio U of the second order temperature fluctuation to the first 
order one for a spherical void (lump) at redshift z, 



U = {AT^^^ /T)ioc{z)/{AT^^yT)iociz), 
which is proportional to c. Plugging Ea. (|5.12p into Eq. (|5.15p . we have 

(C1 + 9C2)' 



U ■ 



W2{y) 



wi{y) 













a \ ^ J 



' (^rn)c 



(5.15) 



(5.16) 
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FIG. 5: The first and second order temperature fluctuations as functions of tlie matter density parameter at present 
rJmO for piiotons passing tfirougfi tiie center of a compensated splierical void at z ~ 0. The matter density contrast at 
the center is (em)c ~ —0.3 and ro = 0.09Hq^, ri — O.lHg^, where Ho is the Hubble constant. The dashed and dashed- 
dotted curves denote (Ar'^-'/r)(oc and (AT''^^ /T)ioc, respectively. The solid curve represents the total temperature 
fluctuation (AT'-'^ /T)ioc + [AT'-^^ /T)ioc. 
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FIG. 6: For a compensated spherical lump. The parameters are as the same as in Fig. (5] The matter density contrast 
is (em)c = 0.3. 
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As shown in Fig. [7l the absolute value \U\ for a fixed value (em)c monotonically increases as the redshift 
z of a void/lump increases. This behavior is naturally expected, since all the models approach to the EdS 

model for which U is infinite (because of (/)'^' = 0) in the limit z ^ oo. It turns out that the ratio U is 
positive (negative) for (em)^ < 0(> 0), respectively. In other words, we have positive (negative) correlation 
for a void (lump) between the first and the second order fluctuations. Thus, we conclude that the first order 
temperature fluctuation corresponding to the linear integrated Sachs-Wolfe (ISW) effect is enhanced (suppressed) 
by the second order effect for a void(lump). 

For (n - ro)/ri = 0.2, or equivalently, y = b/c= 1.049, the behaviors of {AT^^^ /T)ioc{z)/{AT^^^ /T)ioc{0) 
and 0.05 (Ar(2)/r)ioc(^)/(AT(2)/r)ioc(0) are shown in Fig. [51 when c is fixed. In contrast to the first order 
fluctuations which decrease as z increases, it is found that the second order fluctuations have a peak at a certain 
epoch Zp. For a model with {Qmo, ^Ao) = (0.27, 0.73), we have Zp ~ 1. 

VI. CONCLUDING REMARKS 

In this paper we have confirmed that in the nonzero-A cosmological model, the second order solutions (in 
the Poisson gauge) of scalar-type perturbations in general relativistic theory coincide with the solutions in 
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FIG. 7: The ratio between the second and first order temperature fluctuations U for a light path passing through the 
center as a function of the redshift z of the void/lump. The soUd, dashed, and dashed-dotted curves correspond to 
flat models with Qmo = 0.9, 0.6, and 0.3, respectively. The width of the wall is chosen to be 20% of the outer radius, 
i.e., (ri — ro)/ri = 0.2. It turns out that the dependence on the ratio between the width and the outer radius is not 
prominent. 
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FIG. 8: The first order and second order temperature fluctuations as a function of the redshift z of a void/lump with wall 
width (n — ro)/ri = 0.2 , a density contrast parameter c, and density parameters (OmO,f2Ao) = (0.27,0.73). The solid 
and dashed curves denote 0.05 (AT'^^^ /T)ioc{z)/{AT'-^'> /T)io^{0) and {AT'^^^ /T)ioc{z)/{AT'-^^ /T)ioc{0), respectively. 
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the cosmological Newtonian theory in the 'Newtonian' limit, provided that the size of the perturbation is 
sufficiently small compared with the horizon scale. Thus, the second order temperature fluctuation due to a local 
perturbation sufficiently smaller than the horizon scale can be calculated using the Newtonian approximation 
as long as the solutions satisfying the Poisson equation and equations of continuity and motion are used in the 

Newtonian approximation. 

In order to clarify the behavior of second order temperature fluctuations, we considered a simple spherical 
top-hat type local density perturbation (void/lump) in which the pressureless matter is totally compensating 
and whose size is smaller than the horizon size. We have found that the flrst order temperature fluctuation 
corresponding to the linear integrated Sachs- Wolfe (ISW) effect is enhanced (suppressed) by the second order 
effect for a compensated void(lump). As a function of redshift z of the local perturbation, the amplitude of 
the second order temperature fluctiiation due to a void(lump) has a peak at an epoch z ^ 1 for a fixed density 
contrast and the comoving size, whereas the amplitude of the first order fluctuation increases monotonically. 
This means that the second order temperature fluctuations may play an important role especially for quasi- 
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linear objects such as voids with radius r = 100 — 200/i ^Mpc at z '--^ 1, which may lead to some observable 
imprints at an angle ^ 5° in the sky ([l3|, [HI, p7|V 

APPENDIX A: THE MAIN COMPONENTS OF FIRST ORDER AND SECOND ORDER EINSTEIN 

EQUATIONS 



The first order and second order components of Ricci and Einstein tensors are expressed using the symbols S 

1 

and 5asR"^ = {RI)b +6R''^ + 6R''^ and G"^ = R"^ - ^S'^^RTf The energy-momentum tensor is T/; = pu^.u'' . In 
what follows, we assume the metric in (|2.1I) with conditions (|2.8p and (|2.9p . For further details see ([1^, [1^, 

M)- 

For the first order perturbations, we obtain the components of Ricci and Einstein tensors 



a^6Rl = 6 



\ a / a i L a a . 



(Al) 



a'6Rt,^2<^f + 2^^[l\ 



(A2) 



1 



4a" 



-2(- 
a 



'\ 2 



^(1) + A^(i) 



20(1)" + 1^,^(1)' _ A 



and the components of energy-momentum tensor are 



1 z 1 z 



(A3) 



(A4) 



a (5 Jo = -Psa w 
1 



(A5) 



a^iJT^ = 0. 



(A6) 



For the second order perturbations, we obtain the following components 
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'\2 n" 



(^(1))2 



^ = - 4^^(i)^« 4^0(1)^-.^/) - 20(1) + 40(1)^1) + 80(1)0(1)',, 



(A7) 
(AS) 



where 



S Gi = a" 5 G) + Wk, + V'.V + 2(V^'^ - -^^'^ ■ 



0W^«-0«^W+ 30(1)0(1) +40(i)0W^ 



W EE 4 



^^f-2^ 
a ) a 



(^(l))2 _ 8^ 0(1)^(1)' + 0;(1)0(1)' _ 0(1)0(1)' 

1^(1) - 20(/)0(/) - 40(1) A0(i) + 2f0(i) - Vd)] [20(1)" + A^(i) 



0(i)-2V'(i)lV'^' 



(A9) 



(AlO) 
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a^SR^ = 3 



Va/aJ 2L a a . 



r /a'\2 


a" 












a 





(All) 
(A12) 



2 



L2 a la \ a / i 



2a' 
a 



1, 



.(2) , .(2)V 



(2) _ ./,(2) 

1 r 



?0 ^2 



.(2) 



(A13) 



The latter three components S Rq, a S R}), and SGj are those which do not include any second order 
terms consisting of products of first order quantities in S Rq, ci^ SRq, and a? 5G\. For energy- momentum 
tensor, we have 



1 



= Psa 



(A14) 



2 2 



(A15) 
(A16) 



APPENDIX B: SECOND ORDER SOLUTIONS OF EINSTEIN EQUATIONS 



In this Appendix we show the solutions of second order Einstein equations. First, using Eq. (IA9P in Appendix 
A and the expressions of first order perturbations, we obtain 



/iA^(2) + ^^(2)'+r^_f^)'l^(2)+^(2)%^^(2)'_iA 

L2 o LaVo/J a 2 



o 

■ 2 a 



(2)l 



.(2) 



(2)1/ 



2a 



{P'f 



1 / a 
-(1 P 

2 V a 



Xij 

2 



2a' , 



W5,. 



J I] 



a' ,\2 
1 P'\ FF 



I] ■ 



Here FF, 



FjFj and FjFj is divided into following independent terms 



Raj ~I" 



D 



where the following conditions are imposed on Ci and Dij 

S'^ai =0, S'W. 



Functions A, B, C'i and F'i,- are determined as follows 



and 6^''Dij^k 



A T 



AC, 



[200 



F,:AF 



(Bl) 

(B2) 
(B3) 

(B4) 
(B5) 
(B6) 
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Dij — FjFj —"iioSij H — ^0o,ij — (Cij + Cj^i). 

Corresponding to the relation (jB2[) . Eq. (|Bip is divided into four parts 

a a 2 \ a J a 2 



{[(7)" 



2a" 

a 



^{l-^P'y[2A{F^)-FiFi] + 



2 2a' 
a 

2a' 



'\ 2 



7'\ 2 



\ a J L a I \ a J 



a 
a 



J2)' 



2a' 



,(2) 



-2n-^p' 

a 



a 
a 



(2)" 2a^ (2)' _ . (2) 



1 



1 



2 1. a 

Eliminating ^^^^ from Eqs.(|B8| and (|B9l) . we obtain an equation for 0^^^ 

,^9^// 3a' ,/,-,/ r2a" /c 
La V I 

2a' 



" {P' f] 



D, 



r-)'p'+Ui-'-^p' 

\ a J a \ a 

L v a / a a 

(7)t1<-')1>-{ 



- 14 



2 1 r, 2a' 

+ 7 1 

4 I a 

a'\2i , r / 



2«'p' 



--14 - 
a V a 

100 „ _|^_2a^p, 
a 



2 — 



a / a 
/a'N* 
19 - 
a / \ a 

'n2 100 



"\ 2 



(p')2|j7^ 



V a / a J J 



(B7) 



(B8) 
(B9) 
(BIO) 
(Bll) 



(B12) 



Solving this equation and using Eq. (|B9|) . we get the expressions for cf)^"^^ and tp^"^^ given in Eas. (|2.17p - (|2.18p . 

APPENDIX C: THE EXPRESSION OF *„ AND THE INTEGRATION OF [F.rf, *o, AND F 

ALONG A LIGHT PATH 



The term / in Eq.g 
1 



IS 



18 



2 2 
-c r 



for < r < ro 



^rl (c - 36) (c + 6) - [(c + 6)^ {r^/rf - 8(c + 6)6(ro/r)3 - 26^] 



for 7*0 < r < Ti 



100 



*o(0) 



1 
36 



3m2„-4 



cr^ - 6(rf - 4)Yr 



for ri < r, 



where 



^*o(0) 



^r^(c - 36) (c + 6) - lr?[(c + 5)2(ro/ri)« 
L^crl-b{rl-rl)frY'- 



{c + b)b{ro/nf ^2b^ 



For < < tq, we obtain the following integrals by use of u = r/ro and y = b/c 

1 



{F,r)^dx I (-c^r^) = -42/(1 + y) In 
2 



+ 



(«? - ulYl^y\u\ + 2ul) + (1 + y)^.;^{(K? - ulY'Vul - (1 - k^)V2 



(CI) 



(C2) 
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+ u^^ tan" 
1 

■ tan 



.lMK^2)l/2_(i_ ^2)1/2 
2 + 



ut + [{uj - ul){l - ul)] 



; _ "f,;, ^ } + [1 + y(i - - uiy/y^ 



(C3) 



and 



100 
~9~ 



*odx / (IcVg) ^ -42;(1 + y) In + ^^S^^' - |(1 - »;)^/^ [^(1 ~ v'Kl + 2z^^) 
+ (3y - 1)(1 + y)] -liul- ulf'^ \^y\ul + 2ul) - 2y{l + y)/u, - y^u\ 

+ tan-1 ''^2l\7f^'^'2v/''2irrn + + 2^(1 - "?)]'^*"' - 



(C4) 



where ui = ri/rcM* = r*/ro and a; = (r^ — r^)^/^. For r* = 0, we obtain the compensating case (i.e., 
ui^ = l + 1/y) 



f 

Jo 



(Frfdr 



100 
~9~ 



The integration of F is derived as follows : 



^cV3 (l+t/)(l-3yln«i). 



Fdx= I -F^r^-rj) ^^Vdr = c(ro)V(u*) 



(C5) 



(C6) 



where 



(l-M*^) ^ (l+y)(4-u*2^ ' y^Mi^-M*^) +3(l + y)/Mi 



1/2 



1/2 



3(1 + 2/) In 



1/2 



1 + (1 - W,2^ 

In the above derivation we assumed that the mass is totally compensating. For r* = 0, we obtain 

J(0) = -^{l + y) ln(«i) = -^(1 + y) ln(i + 1/y). 



(C7) 



(C8) 
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